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We investigate the effects of rotation on one-dimensional ultracold bosons confined to optical 
ring lattices. First, we show that there exists a critical rotation frequency at which the ground 
state of a weakly-interacting and integer-filled atomic gas is fragmented into a macroscopic 
superposition state with different circulation. Second, we point out several advantages of us- 
ing slightly non-uniform ring lattices. Finally, we demonstrate that different quasi-momentum 
states can be distinguished in time-of-flight absorption imaging and propose to probe corre- 
lations via the many-body oscillations induced by a sudden change in the rotation frequency. 
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1. Introduction 

Macroscopic superposition states, to whicli we will refer as (Schrodinger) cat states, 
are one class of strongly-correlated states which has attracted much interest since 
the early days of quantum theory (Til . Apart from foundational questions ([H they 
have been shown to enable precision measurements at the Heisenberg limit Q and 
serve as resources for various quantum information tasks Q- Efficient cat state 
production is therefore an issue of theoretical and practical importance. 

Among theoretical proposals for creating cat states with Bose-Einstein conden- 
sates one is to load ultracold bosons into rotating ring lattices (fP Hiip . In the 
first part of this paper, we investigate the effects of rotation on the ground state 
properties of bosonic atoms confined in a one-dimensional ring (j3 . For commen- 
surate filling and weak interactions there is a critical rotation frequency at which 
the ground and first excited states are the symmetric and anti-symmetric superpo- 
sition states of opposite circulation. Their energy difference decreases exponentially 
with increasing number of particles . This makes it hard to observe this effect 
even in systems with small numbers of particles. On the other hand, superpositions 
of macroscopically distinct flow have been observed in experiments with supercon- 
ducting quantum interference devices (SQUID) ([Hl-[I7j). This motivates the second 
part of this paper, where we point out that for a slightly non-uniform lattice the 
energy gap scales less severely with the number of particles (.15^. Finally, we discuss 
how one can probe cat-like correlations via many-body oscillations induced by a 
sudden change in the rotation frequency and show that different quasi-momentum 
states can be distinguished in time-of-flight absorption images (\15\f . 

We note that there are related articles on the quantum dynamics of a vortex in 
a Bose gas and topological winding in a ID BEC p7[ ). 
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2. Bosons in rotating ring lattices 

Optical ring lattices can be created experimentally by interfering a Laguerre- 
Gaussian (LG) laser beam with a plane wave co-propagating along the z-direction 
(\18\ [ly^ . By reflecting the combined beam back on itself, one can achieve confine- 
ment along the z-direction, thus creating an stacked array of disk shaped traps. By 
controlling the tunneling between the disks and making it much smaller than the 
corresponding tunneling within each ring, one can implement an array of effective 
ID ring lattices. The lattices can be rotated by introducing a phase shift in one 
the two beams that compose the interference pattern. 

We consider a system of ultracold bosons with mass M confined in a ID ring 
lattice of L sites with lattice constant d. The ring is rotated in its plane (about the 
z-axis) with angular velocity Q. In the rotating frame of the ring the many-body 
Hamiltonian is given by (|77|, 



H 



dx<I> 



'2M 



+ y(x) + 



2M 



(1) 



where a is the s-wave scattering length, V(x.) the lattice potential, the angular 
momentum, and x the 3D spatial coordinate vector. 

Assuming that the lattice is deep enough to restrict tunneling between nearest- 
neighbor sites and the band gap is larger than the rotational energy, the bosonic 
field operator $ can be expanded in Wannier orbitals confined to the first band 
6 = ajT^j(x) with VFj(x) = exp A(x') • dx' W^j(x). Here, VFj(x) are 

the Wannier orbitals of the stationary lattice centered at the site j, A(x) = 0,z x x 
an effective vector potential and aj the bosonic annihilation operator of a particle 
at site j. In terms of these quantities, the many-body Hamiltonian can be written, 
up to on-site diagonal terms which we neglect for simplicity, as ([3 El I^OP 



!)• 



(2) 



Here, hj = a^Sj is the number operator at site j, 9 is the effective phase twist 

, Jj is the hopping energy 



induced by the gauge field, 6 = /^'^^ A(x') • dx' = 
between nearest-neighbour sites j and j + 1: Jj = J dxW^ 
and U the on-site interaction energy: U = ^^jj^ f dx.\Wj 



2M 



V2 + V{x) 



w. 



3. The uniform ring lattice 

We start our discussion with the case of a uniform ring lattice, i.e. Jj = J for all 
j. To understand the effect of rotation on the ground state properties of atoms in 
a rotating ring lattice, it is convenient to write the many-body Hamiltonian ^ in 
terms of the quasi-momentum operators bg = ^ Ylj=i ajC"^'^*''-'/^ where 2'Kq/{dL) 
is the discrete quasi-momentum and q is an integer, g = 0, . . . , L — 1. In this basis 
the Hamiltonian ^ reads 

L-l L-l 

H = -2JH,p + —H,^, = J2Egb% + — blblbi^^g+s-iu (3) 

q=0 q,s,l=0 
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where Eq = —2Jcos[2TTq/L — 9] are single-particle energies and the notation || \\l 
indicates modulo L. The modulus is taken because in collision processes the quasi- 
momentum is conserved up to an integer multiple of the reciprocal lattice vector 
27r/d, i.e. modulo Umklapp processes. 

We begin by considering the non-interacting limit, U = 0. While in a static ring, 
6 = 0, the ground state corresponds to all A'" atoms in the go = quasi-momentum 
state, in a rotating ring the ground state acquires a finite quasi-momentum, '^qo, 
to lower the energy gained by rotation. In other words it becomes a current carrying 
state. If one writes 9 = ^rn + ^ with m an integer and < AO < 2tt, = m 
if < A9 < IT and qq = m + 1 ii ir < A9 < 2tt. The quantum number can also 
be viewed as the winding number, which measures the circulation or "vorticity" of 
the gas: At the critical angular velocity = -^^^xr-j (i-e. A9 = vr), the winding 
number increases by one unit. 

In general, for 17 ^ the ground state is unique and corresponds to a macro- 
scopically occupied state with well defined quasi-momentum: qo = m and qo = m+1 
for 0, < 0,c and 0, > Qc, respectively. However, the fact that exactly at the crit- 
ical frequency, 17 = the single particle energies Eq=rn and Eq=m+i become 
degenerate leads to a (A^ -|- 1)— fold degenerate many-body ground state. The 
A'^ + 1 degenerate energy levels correspond to the states with n atoms m q = m 
and N — n atoms in q = m + 1, with n = 0, . . . ,N. We denote such states as 

In the presence of interactions, one expects the degeneracy to be lifted. For 
weak interactions, i.e. U/J <^ 1, one can use lowest-order perturbation theory 
to account for their effect. In view of the absence of direct coupling between the 
different \n,N — n) states, the energy shifts are diagonal and given by En^ = 
NEm + ^[(A^(A^ — 1) — 2n{n — N)]. Therefore, in lowest order of perturbation 
theory, the states with minimal energy become those corresponding to n = and 
n = N. These two states remain degenerate and higher-order perturbation theory 
is needed to break the degeneracy. At this point one has to distinguish whether 
the number of atoms N is commensurate (CO) or incommensurate (ICO) with the 
number of lattice sites L. 

In the ICO case, there is no coupling between |A^, 0) and |0,A^) and therefore, 
these states remain degenerate for all values of U. In the CO system, on the other 
hand, there are many different paths that couple these two states, and the number 
of paths increases exponentially with N and L. The absence of coupling in the 
ICO case can be intuitively understood by consideringthe total quasi-momentum 
^ — TjW Sg=o^ QbqbqWL- The many-body Hamiltonian (3) exhibits a block diagonal 
form if the quasi-momentum Fock states are ordered according to K. In the CO 
case, = nL with n an integer, the states |A^, 0) and |0,A^) have total quasi- 
momentum K = '^\\mnL\L = and K = ^\\{m + l)nL\L = 0, respectively, and 
thus both of them belong to the K = block. On the other hand, in the ICO case, 
A'' = nL + AA^, and hence the two states |A^, 0) and |0, A'"), belong to different 
blocks as = \\mnL + mAN\L / ||(m -|- \)nL + {m + l)AA^||i. Thus these two 
states are not coupled by interaction. 

In the CO case, we construct an effective 2x2 Hamiltonian by projecting on the 
subspace spanned by the two macroscopic states | A^, 0) and |0, A^) 



H- 



E^'^ A 



2x2 



^0 

A 



E 



(1) 



with A 



n(L-l) 



(4) 
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Where A = Z^,,...^ ^^^^ .^^^^^^^^^Z^y^^E^^ -e^ ^ = ^'^ ™- 
ber density. The Hij are transition matrix elements introduced by the interaction 

term of the Hamiltonian H\at] and Si are either En'^ or non-interacting many-body 
eigenenergies depending upon whether the intermediate states are or are not in the 
N — n) manifold. The factor n{L — 1) corresponds to the minimum number of 
collision processes necessary to connect the states |iV, 0) and |0,A^) and the sum 
is over all the different paths that generate such couplings. At the critical rotation 
frequency we have E^^^ = E^^^ and due to the non-zero value of A the symmetric 
and anti-symmetric superpositions |ib) become the ground and first excited state 
separated from each other by an energy gap 2A 



4. Extension to a ring superlattice 



In the previous section we found that for a uniform ring lattice the coupling A 
between the two macroscopic states scales exponentially with increasing number of 
particles (see Eq. Q), and that the cat state ([s]) only appears in the CO case. 
In this section we explore whether the situation can be improved by introducing 
a lattice modulation, i.e. we generalize our discussion to ring superlattices where 
Jj = J for j even and Jj = t for j odd. 

Due to the superlattice potential the quasi-momentum states \q) and 1^ -|- L/2) 
are coupled and the single-particle Hamiltonian is no longer diagonal in the quasi- 
momentum basis. Via a unitary basis transformation 

/ Cq \ _ ( cos aq 1 siu Uq 

V Cg+L/2 y V i sin cos Uq 



Og+L/2 , 



with 

tan 2aq = ^^-^ • tan ^6* - ^ ) (7) 

we diagonalize the single-particle Hamiltonian i^sp = Y^^=o ^ {^q^q + Eqf^q+L/2) 
where the single-particle energies are given by 



E^ = ±^J-^ + t^ + 2Jtcos(^2e-^y (8) 

The single-particle spectrum Q contains a degeneracy at the critical phase twist 
9 = 11/4 even the presence of the superlattice potential, i.e. we have Eq = E^^, 
but note that for L 7^ 4 these states are not the ground states of the system. 

Following similar arguments as in the previous section we will show that weak 
on-site interactions lift the degeneracy at 6 = 7r/4 and lead to the formation of 
strongly-correlated states in the many-body system. In the absence of interactions 
[/ = the ground state of a bosonic many-body system is the state with all 
bosons occupying the lowest-energy single-particle state. At the critical phase twist 
6 = tt/A the single-particle spectrum has, however, a twofold degeneracy, so that 
there is a A^ -|- 1-dimensional degenerate subspace at A-particle level. A convenient 
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basis for this subspace is given by the Fock states \n,N — n) with < n < A^, where 
n particles are in the single-particle state of energy Eq and N — n particles in the 
one of energy E^^^, respectively. For weak interactions NU/L <^ 2\/ + this 
subspace is the low-energy sector of the many-body problem for all phase twists 9 
and tunneling strength ratios t/J. 
The effective Hamiltonian in first order of the on-site interaction strength U is 

^eff = [EqUo + EljJiLi^ + ^ {2hohL/4 + N'^ - N) + (^-^cIcIcl/4Cl/4 + H.c^ 

(9) 

where the number operators and the parameter rj simplifies for 

e = tt/4 to 7? = (J2 - i^)/(J^ + t^). At the critical phase twist 6 = -k/A the 
terms in the first bracket are an unimportant zero-energy offset, whereas the terms 
in the second bracket shift the energies of the states in the subspace differently, 
e.g. they lead to an energy difference of U{N — 1)/ L between the states |A^, 0) and 
\N — 1, 1), while the states |n, N — n) and | — n, n) remain pairwise degenerate. 
The terms in the third bracket are off-diagonal in the Fock basis of the subspace 
and describe two-particle scattering between the two single-particle modes. They 
lift the remaining pairwise degeneracies in the many-body spectrum. 

Let us now focus on slightly non-uniform rings t/J ~ 1 close to the critical phase 
twist 9 ~ 7r/4. Since the terms in the second bracket of Eq. ^ increase the energy 
for all states in the subspace apart from | A^, 0) and |0, A^) and the coupling between 
the states is weak (as the coupling ry is small in this limit), we project the effective 
Hamiltonian ^ onto the subspace spanned by these two nearly-degenerate lowest- 
energy states (jjjp . As there is no direct coupling between |A^, 0) and |0,A^) we 
calculate the total coupling through intermediate states using perturbation theory. 
After eliminating the intermediate states we obtain, up to a constant energy shift 
which we neglect, the following 2x2 Hamiltonian 

where /S.E is the energy difference between the states |A^, 0) and |0, A^) caused by 
the detuning of the phase twist from resonance IS.9 = 9 — tt/4, i.e. 

and A is the coupling between the states |A^, 0) and |0, A^) due to the off-diagonal 
terms of the effective Hamiltonian As the latter only directly couples the states 
|n, N — n) and |n it 2, A^ — n =p 2), the first non-vanishing order is given by 

A = y>Ol4yiO.^^)^ ^ U (^rlY/' ^^Nl ^^^^ 



where E^^ = {2j{N — j) + N"^ — is the diagonal interaction energy shift. 

In order to understand the influence of the supeiiattice potential on the formation 
of the cat state, we show in Fig. [l|^a), for the example of A'" = 4 particles on L = 4 
lattice sites, the quasi-momentum states \nq=Q,nq=i,nq=2,ng=s) which contribute 
to the ground and first excited state in lowest non-vanishing order of perturbation 
theory. Solid lines signify coupling due to the interaction Hamiltonian i^int which 
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Figure 1. (a) Non-zero matrix elements between the quasi-momentum states \nq—o,nq—i,nq—2,nq—3) 
for N = 4 particles on L = 4 lattice sites. Solid lines signify coupling due to scattering -ffint and dashed 
lines due to the off-diagonal part of the single-particle Hamiltonian Hup in case of the superlattice (J ^ t). 
We only show states which contribute in lowest non-vanishing order of perturbation theory to the cat state 
at the critical phase twist. All states with(out) box have total quasi-momentum K = 2 (K = 0). (b) 
Energy gap as a function of particle number N with t/J = 0.7 (upper curve) and t/J = l (lower curve) 
for the ring of L = 4 sites and U/J = 0.5. The solid lines are exponential fits to guide the eye. 



is present in uniform ring lattices (J = t) as well as ring super lattices (J 7^ t); 
dashed lines stand for the coupling due to the off-diagonal part of the single- 
particle Hamiltonian Hgp induced by the superlattice. The three states |3,0, 1,0), 
|0, 3, 0, 1) and |2, 2, 0, 0) have total quasi-momentum K = 2 and are not coupled in 
the case of the uniform ring lattice. They introduce four efficient coupling paths 
involving only N/2 = 2 scattering processes and the intermediate state |2,2,0,0) 
lies in the degenerate subspace. 

Let us compare our results for the coupling A in the case of a uniform ring lattice 
^ with the case of a ring superlattice (Il2|) . In both cases the coupling A decreases 
exponentially with the number of particles N . This is because multiple two-particle 
scattering processes are the microscopic origin of the coupling A and transitions 
between the two configurations |iV, 0) and |0, A/") are thus highly off-resonant N/2- 
th and n{L — l)-th order processes, respectively. Consequently, the energy gap 
vanishes in the thermodynamic limit — > 00 and cat state production both in 
uniform and slightly non-uniform ring lattices is restricted to modest numbers of 
atoms. However, as the energy gap in first non- vanishing order is proportional to 
C/(r;/2)^/2 foj. slightly non-uniform rings as compared to C/(C//J)"^'^ for uniform 
rings, it can be one order of magnitude bigger in the relevant parameter regime. 
Moreover, the appearance of cat states in slightly non-uniform rings is not limited 
to commensurate filling as long as the number of particles N is even. 

In Fig. [ij^b) we plot the energy gap L = 4 and U/J = 0.5 as a function of number 
of particles N for t/J = 0.7 and t/J = 1. The solid line is an exponential fit to 
the points to guide the eye. We see that by changing from t/J = 1 to t/J = 0.7 
the energy gap increases by over one order of magnitude for numbers of particles 
> 8. As experiments are typically limited to 1 sec and typical hopping energies 
are of order of 0.05Er ~ lO^Hz, the detection of cat states by non-equilibrium 
dynamics is limited to A/ J = 10~^. This in turn restricts our cat-state production 
scheme to a few tens of particles in contrast to less than ten particles for uniform 
ring lattices (see Fig. [T|^b)). 



5. Dynamical detection of cat-like correlations 

In this section we propose to induce many-body oscillations by suddenly changing 
the applied phase twist 9 in order to detect the coherent superposition of two quasi- 
momentum states |A^,0) = (bl)^ \vac) / and |0, A^) = {b{^^)^\vac) /Vn\ at the 
anti-crossing of the many-body spectrum. 
If we assume the system is initially in the ground state of the full Hamiltonian at 
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Figure 2. (a) Many-body dynamics |6o(i)P and |f)jv(*)P for L = 4 and N = 8 following a sudden change 
in phase twist from 9 = to 9 = 7r/4 at U/J = 1 and t/J = 0.7. (b) Transverse time-of-flight images of 
an initial state with quasi-momentum ijo = and qo = L/A for L = 16. 



= which will predominantly be the state [A'', 0), i.e. \'ip{t = 0)) « [A'', 0). Then, 
the phase twist is changed to 9 = 7r/4, where the eigenstates are approximately 
|ib) (|A^, 0) ± |0, N))/^/2. The time evolution of the probability to be in the state 
\N,0) and \0,N) is 

Mtr = \{N,0\m)\'-''-^^ (13) 

and |6Ar(t)p 1 - |6o(i)P where hu = y^AE'^ + |2A|2 denotes the energy gap and 
where the approximate signs indicate that we neglect the depletion due to the non- 
diagonal part of the single-particle Hamiltonian as well as the on-site interaction 
Hamiltonian. 

In Fig. |2|a) we show the many-body dynamics, i.e. |6o(i)P and \bN{t)\^ for L = 4 
and N = 8 following a sudden change in phase twist from = to ^ = 7r/4 at 
time t = with U/J = 1 and t/J = 0.7. We see that the many-body oscillations 
following this sudden change occur at the frequency of the gap hv. These oscillations 
are modulated by an oscillation with frequency U{N — 1)/L which is the energy 
difference between the nearly-degenerate ground states and the next lowest-lying 
intermediate states \2,N — 2) and — 2,2). The amplitude of the oscillations is 
smaller than one, i.e. |6o(i)P + I^Af(*)P < due to finite depletion of the cat state. 

In Fig. we show the numerically calculated time-of-flight images for two 
states with initial quasi-momentum qo = and L/4, respectively. These states 
correspond to the two types of initial distributions that one has to experimentally 
distinguish in order to probe the many-body oscillations. The figure shows the 
distinctive interference pattern of the two initial states. 

In the limit of L ^ oo lattice sites the time-of-flight profile reduces to 

|V'(r,i oo)p oc 

where Jqa{x) are Bessel functions of the first kind and R denotes the radius of the 
ring lattice. This case corresponds to the rotationally symmetric case discussed in 
Ref. I\21\\ in which a state initially with zero quasi-momentum for example will ex- 
hibit an interference peak at the origin while a state with non-zero quasi-momentum 
will exhibit a central hole. For a finite number of lattice sites the sum does not 
correspond exactly to a Bessel function and the momentum distribution does not 
become fully radially symmetric. However, there is still a unique correspondence 
between the position of the peaks in the absorption image and the initial quasi- 
momentum distribution (see also Ref. i^))- Consequently, the latter provides a 
means to experimentally determine the quasi-momentum of the wavefunction be- 
fore the release. 



J 



Mr 
htR 



(14) 
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6. Conclusion 



We have investigated the effects of rotation on the ground state properties of 
bosonic atoms in a rotating ring lattice. In the uniform geometry only the com- 
mensurate (CO) system exhibites a cat state at the critical phase twist and its 
gap scales with the number of particles N as {U/2L)"^^~^\ Introducing a super- 
lattice lifts the CO constraint and improves the scaling to C/(?7/2)^/^. This makes 
the cat less sensitive to imperfections in the system. Finally, we have proposed to 
probe the cat-like correlations via the many-body oscillations induced by a sudden 
change in the rotation frequency. Since the exponential scaling of the energy gap 
with the number of particles remains, we suspect that cat state production in ul- 
tracold atomic gases is limited to modest number of atoms at least for systems with 
contact interactions. One possible way-out is to consider systems with long-range 
interactions such as polar molecules with dipolar interactions. 
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